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The influence of an asymmetric in-plane magnetic anisotropy K x ^ K y on the thermally acti¬ 
vated spin current is studied theoretically for two different systems; (i) the F\N system consisting 
of a ferromagnetic insulator (F) in a direct contact with a nonmagnetic metal (IV), and (ii) the 
sandwich structure iV| JP’IJV consisting of a ferromagnetic insulating part sandwiched between two 
nonmagnetic metals. It is shown that when the difference between the temperatures of the two non¬ 
magnetic metals in a A5|j^lAf structure is not large, the spin pumping currents from the magnetic 
part to the nonmagnetic ones are equal in amplitude and have opposite directions, so only the spin 
torque current contributes to the total spin current. The spin current flows then from the nonmag¬ 
netic metal with the higher temperature to the nonmagnetic metal having a lower temperature. Its 
amplitude varies linearly with the difference in temperatures. In addition, we have found that if the 
magnetic anisotropy is in the layer plane, then the spin current increases with the magnon temper¬ 
ature, while in the case of an out-of-plane magnetic anisotropy the spin current decreases when the 
magnon temperature enhances. Enlarging the difference between the temperatures of the nonmag¬ 
netic metals, the linear response becomes important, as confirmed by analytical expressions inferred 
from the Fokker-Planck approach and by the results obtained upon a full numerical integration of 
the stochastic Landau-Lifshitz-Gilbert equation. 

PACS numbers: 

I. INTRODUCTION 

One of the key observations that gave impetus to the 
field of spin caloritronics was the discovery of the spin 
Seebeck effect (SSE)i which amounts to the emergence 
of a spin current upon an externally applied thermal 
gradient^— . Technologically, various SSE-based nano¬ 
electronics devices are envisaged. For example, portable 
thermal diodes have been proposed to control and rectify 
the heat and spin currents. 

The SSE was observed in materials of substantially 
different transport properties such as metallic ferromag- 
net Co 2 MnSi, semiconducting ferromagnet GaMnAs, and 
magnetic insulators LaY 2 Fe 50 i 2 and (Mn, Ze)Fe 2 C> 4 . 

Thus, the underlying mechanism may well depend on the 
specific case under study as well as on the experimen¬ 
tal setup. In metallic ferromagnetic systems, the spin is 
transferred via charge carriers activated by the thermal 
bias, while in magnetic insulators the SSE is mediated 
by magnons flowing towards the cold edge of the sam¬ 
ple. A theory of the magnon driven SSE was developed 
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in Ref. [flj and implemented beyond the linear response 
regime Thereby, the concept of magnon tempera¬ 
ture is of a key importance for understanding the physi¬ 
cal origin of the spin current in magnetic insulators: An 
external heat bias applied to the system thermalizes the 
phonon subsystem much faster than the magnons relax. 
Therefore, the magnon temperature T™ is influenced by 
the already established phonon temperature profile Tp. 
The thermally activated spin current is related to the 
temperature difference between the phonon and magnon 
subsystems. 

Recent studies based on the macrospin approach (valid 
for samples of small dimensions on the range of the ex¬ 
change length) concern the linear—, and nonlinear re¬ 
sponse regimes^. In both regimes, the obtained ana¬ 
lytical expressions for the spin current are proportional 
to the difference between the phonon and magnon tem¬ 
peratures. As shown in Ref. [LJ, the result for the spin 
current obtained in the linear response theory is a par¬ 
ticular case of the result obtained in the Fokker-Planck 
approach, and corresponds to the low magnon tempera¬ 
ture regime. The nonlinear effects substantially change 
the role of the magnetic anisotropy in the formation of 
the spin current. In the linearized approach, the role of 
the magnetic anisotropy is similar to that of an external 
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magnetic field, and can be described by a certain effective 
field. This is not the case when magnetic fluctuations are 
large which is more likely for higher temperatures. 

A quantitative criterion for the threshold magnon tem¬ 
perature, above which the anisotropy plays an impor¬ 
tant role, is defined by the following inequality^, T™ > 
M s VH e ff/kB- Here M s is the saturation magnetization, 
V is the total volume of the ferromagnet, and H e s is the 
effective magnetic field. In the present work we focus 
on phenomena inherent to the nonlinear regime. In par¬ 
ticular, we study the influence of the in-plane magnetic 
anisotropy on the SSE. We show that the effect of the in¬ 
plane anisotropy on the thermally activated spin current 
is different from the effect of the uniaxial out-of-plane 
anisotropy studied in Ref. 0 

We consider two different system relevant for the lon¬ 
gitudinal SSE: (i) an F\N structure consisting of a fer¬ 
romagnetic insulating part (F) attached to a nonmag¬ 
netic metallic part (N), and (ii) a ferromagnetic insulat¬ 
ing part sandwiched between two nonmagnetic metallic 
ones (iV|T|./V). In the following the nonmagnetic metal¬ 
lic part will be referred to as metallic part or simply as 
metal, while the ferromagnetic insulating part will be 
referred to as ferromagnetic part or simply as ferromag¬ 
net. We show that in the both cases thermally activated 
spin current is parallel to the temperature gradient. In 
the case of N\F\N structure we show that if the differ¬ 
ence between the temperatures of the two metals is not 
large, and the temperature dependance of the spin con¬ 
ductance can be ignored, then the spin pumping currents 
from the ferromagnet to the adjacent metals are equal in 
amplitude and are oriented oppositely, so they do not 
contribute to the total spin current. The only contri¬ 
bution to the total spin current is then due to the spin 
torque current. We show that the spin torque current 
flowing through the ferromagnetic part is a linear func¬ 
tion of the difference between the temperatures of the two 
metals. Though, the considerations apply directly to the 
case when the ferromagnetic part is insulating, the re¬ 
sults are also applicable when the ferromagnetic part is 
metallic. 

In Section 2 we present the model F/N structure. The 
Fokker-Planck approach is briefly described in Section 3, 
where analytical results for the spin current in the F\N 
system are presented. Numerical results for the total spin 
current in the F\N system, obtained by a direct numeri¬ 
cal integration of the stochastic Landau-Lifshitz-Gilbert 
(LLG) equation, are presented and discussed in Section 
4. We show that in the presence of an in-plane magnetic 
anisotropy and a weak magnetic field, the (F\N) system 
supports a spin pumping current only if the axial sym¬ 
metry in the system is broken by an in-plane magnetic 
anisotropy. In Section 5 we discuss analytical results 
on the spin current flowing through the system N\F\N, 
while the corresponding results obtained by numerical 
integration are presented in Section 6. We show that 
such a system supports spin pumping current only for a 
sufficiently large thermal gradient. Summary and final 


conclusions are in Section 7. 


II. MODEL OF THE F\N STRUCTURE 

We consider first the thermally activated spin current 
through the interface in the F\N system. Our main focus 
is on the influence of the in-plane magnetic anisotropy, 
-M s ( K x ml + Kyrriy)/ 2, on the spin current. Here, m is 
a unit vector along the magnetic moment, m = M/M s . 
We assume that the thermal equilibrium between the 
electron and the phonon subsystems in the metallic part 
as well as in the ferromagnetic part is restored internally 
much faster than the equilibrium between the two parts. 
In terms of the local temperature, which is based on the 
hierarchy of relaxation times, this means that the tem¬ 
peratures of the phonons, I'n(p) > an< ^ the electrons Tn( F ) > 
baths are equal in both the metallic and ferromagnetic 
parts, Tjy = Tp = T/v, Tp = Tp = Tp. However, there 
is a difference in temperatures of the two components, 
Tp p T/v, which can be controlled externally. This dif¬ 
ference drives the SSE. The interaction between the non¬ 
magnetic and ferromagnetic subsystems is mediated via 
the magnon bath, described by the magnon temperature 
Tp. Due to the slower magnon relaxation, Tp may be 
different from Tp. 



FIG. 1: A schematic illustration of the considered F\N sys¬ 
tem. The system consists of a ferromagnetic insulator (left) 
and a nonmagnetic metal (right) in direct contact. The tem¬ 
peratures of the ferromagnetic (Tp) and metallic (T/v) parts 
are in general different. 


The magnetization dynamics of the ferromagnetic part, 
to be considered in the following section, is described 
by the LLG equation in the macrospin approximation. 
Such an approach excludes nonuniform magnetization, 
and therefore is applicable when the magnetic compo¬ 
nent of the system has a small volume, i.e., its lateral 
and vertical dimensions are small, usually in the nanome¬ 
ter range. Using the Fokker-Plank-equation technique, 
we evaluate the mean value of the spin current flowing 
through the interface from the ferromagnetic part to the 
metallic one. We also study there the dependence of the 
total spin current on the in-plane magnetic anisotropy. 
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III. FOKKER-PLANCK APPROACH: 

ANALYTICAL SOLUTION FOR THE SPIN 
CURRENT IN THE F\N SYSTEM 

The total spin current flowing through the interface 
in the F\N system consists of two contributions - the 
spin pumping current flowing from the ferromagnet to 
the normal metal, and the thermally activated spin cur¬ 
rent flowing in the opposite direction (in the following 
also referred to as the spin torque current). Powered 
by the phonon bath, thermal noise leads to formation of 
the fluctuating spin torque current in the normal metal, 
I flit). Effect of this fluctuating spin torque current on 
the ferromagnetic insulator can be described by a ran¬ 
dom magnetic field h'{t) acting on the magnetization^ 
Iji{t) = — I ' 1 ^ 7 m{t) x h'(t). Here M s is the satura¬ 
tion magnetization, V is the total volume of the fer¬ 
romagnet, and 7 is the gyro-magnetic factor. On the 
other hand, the thermally activated magnetization dy¬ 
namics in the ferromagnet gives rise to a spin pumping 
current emitted from the ferromagnet into the normal 
metal, I s = ^f^-a'fn x to. Thus, the total average dc 
spin current across the interface can be written in the 

form&M 


(Is) = 


M k V r 


7 


a (to x to) — t(to x h!) 


(1) 


where , a' = (' yh/ATrM s V)g r is the magnetization damp¬ 
ing constant (related to the spin pumping), and g r is the 
real part of the dimensionless spin mixing conductance. 
Furthermore, h'{t) is the random magnetic field. If the 
random thermal force has correlation time much shorter 
than the response time of the system, one can assume 
that the noise is white. A quantitative criterion for using 
a white noise is£ k B T hio o, where loq is the ferro¬ 

magnetic resonance frequency. Taking into account the 
fact that the response time of the system is defined by 
uq ~ 1 GHz, for the low temperature limit one obtains 
T >> 0.01K. Evidently, this temperature regime allows 
using the white noise in our proble m 19 ’ 20 : 

( 7 h'(t) 7 h'(t'))=a'%6(t-t'), ( 2 ) 

for i,j = x,y,z, where a' 2 = 2a! 7 ksT^/M S V. Thus, to 
find the spin current we need to know the time evolu¬ 
tion of the magnetization. The magnetization dynamics 
is described by the stochastic LLG equation for the di¬ 
mensionless unit vector m 

to = - 7 m x I H e ff + h ) + am x to. (3) 


The effective magnetic field U e ff consists of the external 
field applied along the z axis and of the in-plane mag¬ 
netic anisotropy field Ha = —<914 n/dM = K x m x i x + 
K y m y i y , where i x / y ) are the unit vectors along the axis 
x(y). The magnetic anisotropy energy density, U an = 
— \M s {K x m 2 + K y m 2 ), is described by the anisotropy 


constants K x and K y . In the above equation, a is the 
total magnetic damping constant, aT™ = a^Tp + a'Tp[. 
This damping constant includes the contributions from 
the standard bulk damping constant op which is associ¬ 
ated with the lattice thermal oscillations (Gilbert damp¬ 
ing) and from the damping constant a! associated with 
the contact to the normal metal. Introducing the en¬ 
hanced total damping constant a has a clear physical mo¬ 
tivation. Due to the effect of the F\N interface, magne¬ 
tization dynamics in the F layer is additionally damped, 
and this enhanced damping is due to a magnonic spin cur¬ 
rent transferred from the ferromagnetic insulator to the 
normal metaL 21 Also here we assumed that the random 
contributions from the uncorrelated noise sources are 
totally independent and, therefore, the total enhanced 
damping constant is factorized^. Finally, h in Eq.(3) is 
the total random field, with the corresponding correla¬ 
tion function of the form^ 

(. hi(t)hj(t ')) = a 2 dijS{t - t'), (4) 

where a 2 = 2akBT™/M s V 7 is the coefficient propor¬ 
tional to the magnon temperature. 

Following the procedure described in Ref. 13 we de¬ 
rive the Fokker-Planck (FP) equation for the distribution 
function /(to, T): 


di _ 1 <9 
dt 1 + a 2 dm, 


(rh x Weft) / 


+ am x (to x w e fi) / — 


2(1 +a 2 ) 


TO X TO X 


df_ 

dm, 


( 5 ) 


where w e ff = 7 H e g = {u> x m x ,ujymy,uJo). The stationary 
solution of the FP equation reads: 

/(to) = Z~ 1 exp (b J ujes-drhj , ( 6 ) 

where Z = J exp (/3 J w e ff ■ dm) d 3 m is the normalization 
factor, and we introduced the following notation: /3 = 
2a (l + a 2 ) /a 2 . Taking into account Eqs. dTJ to (0, 
after some cumbersome, but straightforward calculations 
we find the following expression for the average total spin 
current: 


(Isz) = a'k B [T^A{ 1 - TO 2 ) - 2Tp B x (m z m 2 ) 
(m z m 2 ) - 2 Tjv<to z >}, 


o rpm TD 

p JDy 


( 7 ) 


where A = B r = 


_ arif] o _ KyMsV 

KB J- F 2k B Tp l > D y 2k B T™ ’ 

while the averages occurring in Eq. Q are defined in the 
Appendix [see Eqs. (1A1I) to (1A5[) ]. By definition, the spin 
current is generally a second rank tensor and is charac¬ 
terized by spatial orientation and projection of momen¬ 
tum, I x ^y ’z’ z ■ However, due to the particular geometry 
under consideration, only the longitudinal spin current 
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component, I^ z = I sz . is nonzero, while the transversal 
spin current components vanish, I sx = I sy = 0. Thus 
our findings are in favor of the recent experiment^ in 
which upper limit for the transverse spin Seebeck ef¬ 
fect was observed several orders of magnitude smaller 
than previously reported. In spite of the fact that ex¬ 
pression of the spin current doesn’t depend on the en¬ 
hanced damping constant a explicitly, due to the rela¬ 
tion aT™ = a 0 T F + a'T N spin current still depends on a 
implicitly. One can invert dependence on the spin mix¬ 
ing conductance damping a! into the dependence on a. 
However precise values of a' is too much related to the 
characteristics of particular interface. Therefore for the 
sake of general interest we quantify spin current in terms 
of < I sz > /oiks- The in-plane magnetic anisotropy has 
different physical consequences, when compared to the 
case of an out-of-plane anisotropy studied in Ref. [3 
First of all, the expression for the total spin current in 
the case of an in-plane anisotropy is different from that 
for an out-of-plane anisotropy. Only in the case of an 
axial symmetry, K x = K y , the expressions for the to¬ 
tal spin current become partly similar. However, as we 
will see below, the main effect of the in-plain anisotropy 
concerns the asymmetric case, K x ^ K y . 

The expression for the total spin current, Eq. Q, is 
quite general. Therefore, we consider now some asymp¬ 
totic situations. In the symmetric case, K x = K y = Ha, 
one finds B x = B y = B = HaM s V/2k B Tp . Equa¬ 
tion (ED reduces then to a simpler form, while the mean 
components of the magnetization can be calculated ana¬ 
lytically: 


/ i A 

(m »> = 2 B 


e~^~ B 


sinh[A] 


\J 7 tB G(A,B) 


( 1 -to 2 ) = 1 - — - — 


2 B 


2 B 


e 4 ^ b (2B cosh(A) + 74sinh(A)) 


(m z (r 


G(A,B) 

v. __/_AV_3_A 

~ 1 OR ) 4 52 ^ 1 R 


2 B J 


A R 

e ~TW~ B 


(8a) 


(8b) 


2 B 


2 v / 7rH 5 / 2 G(H, B) 

x [2ARcosh(A) + (A 2 + 4 B) sinh(A)] . (8c) 

Here, we introduced the following notation: G{A,B) = 



erf(- • •) is the error function. Interestingly, the in-plane 
magnetic anisotropy in the symmetric case is equivalent 
to the out-of-plane anisotropy, — Haiti 2 . Therefore, after 
substituting B —>■ — B in Eqs ©, we recover the results 
derived earlier—. 

In the case of weak magnetic anisotropy, B < 1, the 
mean components of the magnetic moment can be further 
simplified (see Eq. (IA6I) to Eq. (IA9I) in the Appendix). 
In turn, in the absence of magnetic anisotropy, B —> 



FIG. 2: Dependence of the average total spin current 
{Isz}/ot ks on the external magnetic field Ho for the fol¬ 
lowing parameters: K x = K y = Ha > 0 (with Ha as in¬ 
dicated), T 'P = 300 K, Tjv = 297 K, M s = 800 G, and 
V = 1.6 x 10 -18 cm 3 . The figure shows that the spin cur¬ 
rent decreases with the increasing the magnetic anisotropy. 
Since we measure spin current in the units of {Isz}/ot ks spin 
current doesn’t depend on the a'. Besides analytical result 
is obtained via the steady sate solution of the FP equation. 
This implicates that system is already relaxed. This is reason 
why spin current doesn’t depend on the a as well. Values of 
the inverse damping constant 1 /a defines relaxation rate and 
steady sate distribution function is formed beyond this time 
scale. 


0, from Eqs. ED and © one finds (m 2 ) = L{A), ((1 — 
ml)} = 2 L(A)/A, and (J S2 ) = 2a'k B L{A){T™ - T N ), 
where L(A) = coth(A) — 1/A is the Langevin function. 
This result (obtained for B —> 0) recovers the previously 
obtained expression for the total spin currenl— . 

Now we address the case of a high magnon tem¬ 
perature, M s VH 0 /k B T™ < 1 and H A M s V/2k B T™ < 
1. Then, from Eqs. ED and © we obtain {I sz ) = 

2 a 'k B ^k^ ( T f ~ T N ) (l - fOT the total 

spin current. From this formula follows that the posi¬ 
tive in-plane anisotropy {Ha > 0 ) reduces the total spin 
current, while the negative in-plane anisotropy {Ha < 0 ) 
enhances the total spin current. Apart from this, the spin 
current is proportional to the difference (T™ — Tn ) be¬ 
tween the magnon temperature and the temperature of 
the metal. The amplitude of the total spin current, how¬ 
ever, is reduced because of the anisotropy term 
Another interesting observation is that the spin current 
vanishes, (J sz ) —> 0, for zero magnetic field, H 0 > 0. 
This result is rather clear since the dynamics of the mag¬ 
netization is then strongly exposed to the thermal fluc¬ 
tuations, and therefore suppresses the spin current gen¬ 
eration. 


In the case of a strong magnetic field and a weak 
anisotropy, \HaM s V < k B T ™ < H 0 M S V , we enter a 
different regime, and the expression for the total spin 
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FIG. 3: The average total spin current as function of the 
external magnetic field Ho for the following parameters: K x = 
K y = Ha < 0 (with Ha as indicated), Tp = 300 K, T/v = 
297 K, M s = 800 G, and V = 1.6 x 10 -18 cm 3 . Absolute 
values of the magnetic anisotropy field are the same as in 
Fig. [2j but here they have the opposite sign. Note, the spin 
current now increases with increasing \Ha\- 


current in this regime reads 


(Isz) 


2 a'k B (TF 



3 H a k B T ™ \ 
H 0 M S VH 0 ) 


(9) 


The first term in Eq. 0 is the standard contribution to 
the spin current in the linear response, while the second 
term is due to the magnetic anisotropy. As before, the 
positive in-plane anisotropy (Ha > 0 ) suppresses the to¬ 
tal spin current, while the negative in-plane anisotropy 
(Ha < 0 ) enhances the spin current. 


Having considered some limiting asymptotics for the 
thermally activated spin current, let us go back to the 
general solution, see Eqs. 0 and ®. Since the gen¬ 
eral solution is relatively complex for its illustration we 
plot the total average spin current as a function of the 
external magnetic field Hq and the magnon tempera¬ 
ture T™. First, we consider a symmetric situation, 
K x = K y = Ha, and then proceed with more impor¬ 
tant an asymmetric case, K x ^ K y . 

Figures [2] and [ 3 ] show the dependence of the total spin 
current on the external magnetic field Hq applied along 
the z axis. As one can note, the magnitude of spin cur¬ 
rent depends on the sign of the anisotropy field Ha- For 
positive magnetic anisotropy, Fig. [21 the spin current de¬ 
creases with increasing the anisotropy field, while for neg¬ 
ative anisotropy, the spin current grows with increasing 
absolute value of Ha, see Fig. [3] In turn, the depen¬ 
dence of the average spin current on the magnon tem¬ 
perature is presented in FigUfor the indicated values of 
the anisotropy constant K x and the constant K y . All 
curves cross at the point T™ = Tn, where the spin cur¬ 
rent is equal to zero as the system is then in thermal equi¬ 
librium. Interestingly, the spin current for K x > K y is 


FIG. 4: The average total spin current as the magnon tem¬ 
perature T™ varies. The following parameters are chosen: 
K y = 100 Oe, H 0 = 200 Oe, T N = 300 K, M s = 800 G, 
V = 1.6 x 1CP 18 cm 3 , and K x as indicated. 


reduced in comparison to the spin current in the symmet¬ 
rical case (K x = K y ). This is because the x-component 
of magnetization increases with the increasing constant 
K x (positive), and thus the magnetization projection on 
the z-axis becomes reduced. 


Before proceeding to a more important case of asym¬ 
metric in-plane anisotropy (K x ^ K y ), we plot the de¬ 
pendence of the mean magnetization component (m z ) on 
the applied magnetic field. As follows from Fig. 0 the 
magnetization component (jriz) increases with the ap¬ 
plied magnetic field, approaching the maximum (m z ) « 
1. However, the limit (m 2 ) = 1 corresponds to the zero 
magnon temperature, A = —> oo, T™ = 0, and 

therefore is beyond the Fokkcr-Plank approach. 


Finally, let us consider the limit of low magnon temper¬ 
ature. Taking into account Eq. 0, Eq. 0, and applying 
the saddle-point method, we find the following formula 
for the average total spin current in the limit of a weak 
magnetic field 2 H 0 / (K x + K y ) < 1: 

Here, we introduced the following notation: r]i = 
H 0 / (I\ x + K y ) and ry 2 = (K x - K y ) (l - 4 rfi). In 
the opposite case of 2H 0 / (K x + K y ) > 1, the spin cur¬ 
rent is given by the expression ( I sz ) = — 2 a'fcsT/v. 

The second term in Eq. m corresponds to the 
spin torque current flowing from the metal to the 
ferromagnet, so it is negative, —oi'IcbTIiTn- Taking 
into account the parity of Bessel functions Ii^/T^ 1 ) 
and I 0 (rj 2 /T™), it is easy to see that the first term 
—a'ksTiiV 2 I 1 (V 2 /T™)/Io(r] 2 /T™) is even with respect to 









6 



FIG. 5: Dependence of the averaged component of magne¬ 
tization \Tn z ) on the external magnetic held Ho for the fol¬ 
lowing parameters: K x = K y = Ha > 0 (with Ha as in¬ 
dicated), Tp = 300 K, Tv - 297 K, M a = 800 G, and 
V = 1.6 x 10 -18 cm 3 . The magnetic held Ho = (0,0, Ho) 
tends to align magnetization vector along the Z axis, while 
the in-plane magnetic anisotropy favors the in-plane align¬ 
ment of the magnetization. 

the permutation (K x — K y ) — > — (K x — K y ), and is al¬ 
ways negative. It disappears in the symmetric case of 
K x = K y while in the antisymmetric case it decreases 
the spin torque current, see Fig. [4] in the low magnon 
temperature regime. 

IV. NUMERICAL SOLUTIONS FOR THE F\N 
SYSTEM 

The results derived analytically in the previous section 
from the solution of the Fokker-Planck equation can be 
further confirmed by exact numerical integration of the 
stochastic LLG equation, assuming the random field h 
as a Gaussian white noise defined through the correlation 
function, see Eq. ©• The magnon temperature Tp is im¬ 
plemented into the simulations via the strength parame¬ 
ter of the correlation function, a 1 = 2aksT™/MgVj (see 
also section 3). 

To solve the stochastic LLG equation we used the Heun 
method, which converges in the quadratic mean to the so¬ 
lution interpreted in the sense of Stratonovich^. From 
the solutions of the stochastic LLG equation we gener¬ 
ated the random trajectories for a sufficiently long time 
interval, until the magnetization components reached the 
stationary regime. This procedure has been implemented 
many times in order to construct an ensemble of the ran¬ 
dom solutions of the stochastic LLG equation^. Each of 
these random solutions corresponds to a certain realiza¬ 
tion of the random noise, while the statistical average 
over the ensemble of realizations designates the mean 
values of the magnetization components. Such average 
components of the magnetization vector (m and m) can 
be used afterwards for the evaluation of spin current. 
This numerical procedure is in general computationally 



FIG. 6 : The total spin current as a function of the exter¬ 
nal magnetic field Ho, obtained via numerical integration 
of the stochastic LLG equation for the indicated values of 
the positive anisotropy constant Ha and for the parameters: 
TP = 300 K, T n = 297 K, M s = 800 G, V = 1.6 x 10 " 18 cm 3 , 
a = 1.0 and a' = 0.05. 

expensive^, since the number of realizations needed to 
reach a good accuracy of the solution to the stochastic 
LLG equation is about one thousand^. 



FIG. 7: Total spin current as a function of an external mag¬ 
netic field Ho, obtained via a numerical integration of the 
stochastic LLG equation for the indicated values of the neg¬ 
ative anisotropy constant, and for the parameters: Tp = 
300 K, Tjv = 297 K, M s = 800 G, V = 1.6 x 10 ” 18 cm 3 , 
a = 1.0 and a' — 0.05. 


This numerical procedure can be used to calculate the 
spin pumping current from the ferromagnetic part to the 
metallic one (first term in Eq. dTJ). In order to calculate 
the spin torque current from the metal to the ferromag- 
net (second term in Eq. dTJ), we have to consider the 
random magnetic field in the metallic part, h', and the 
corresponding correlation function, Eqs. ©. The metal 
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temperature T'n is interlaced with the strength of this 
stochastic field, 2a'k B T N /M S V^. 

Results of the numerical simulations are presented in 
Figs. [ 6 ] and 0 As we see, the numerical results for the 
total spin current are very close to the results obtained 
by means of the Fokker-Planck equation, see Fig. [2] and 
Fig. [3] In both cases, the total spin current increases 
with an external magnetic held. In the case of a pos¬ 
itive anisotropy Held, the larger is the anisotropy, the 
smaller is the spin current, while in the case of a negative 
anisotropy Held, the spin current increases with increas¬ 
ing the anisotropy field. 


V. SPIN CURRENT IN N\F\N STRUCTURES 

The technique used above can be also employed to cal¬ 
culate the thermally-induced spin current in a A|F|iV 
system, shown schematically in Fig. [5] Now, the total 
spin current includes four terms, 

Is — Ifl 1 T Ispl Ifl2 T Isp2i (1-Q 

where I sp i and I sp 2 stand for the spin pumping current 
from the ferromagnet to the metallic parts N 1 (left) and 
N2 (right), respectively (see Fig. [ 8 ]). In turn, Ifn and Ifi 2 
stand for the spin torque current flowing, respectively, 
from the left and right metallic parts to the ferromagnetic 
one due to thermal fluctuations. We assume that the two 
metals have generally different temperatures, Xjvi and 
Tn2- 



FIG. 8: Schematic presentation of the AIEIIV structure con¬ 
sisting of a magnetic element sandwiched between two non¬ 
magnetic parts - N1 on the left and N2 on the right side. 
Temperatures of the parts IV1 and N 2 are generally different. 


Upon laborious calculations, one finds the components 
of the spin pumping current flowing from the ferromagnet 
towards the two (* = 1 , 2 ) metallic parts, 

( I$)=-—a'(r Nl )f{T™), 

M V (12) 

< J ip>) = -^a'(T N2 )f(T?), 

where, for shortness, we introduced the following func¬ 
tion: 


In turn, for the spin torque components we find 
(I^) = 2k B {m z )a'{T N1 )T N1 , 

(14) 

(iff) = -2k B {m z )a'(T N2 )T N2 . 

As one can see from Eqs. (13) and (14), the differ¬ 
ence in the two components of the spin pumping cur¬ 
rent transferred from the ferromagnetic part to the met¬ 
als, (/j^) and (. IsV), is related to the temperature de¬ 
pendence of the damping constants, or, more precisely, 
to the temperature dependence of the spin conductance, 
a'(Tpj ) = 47r Xf v (Tn) • Such a temperature depen¬ 
dence of the spin conductance g r (TN) has been measured 
in a recent experiment^*. 

For convenience, we denote c/(Xjv 1 ) = a! and 

a'(Tff 2 ) = ex' + Ac/, and rewrite the expression for the 
total spin current in the form 


(I sz ) = 2 k B {m z (T™))a'(T N1 - T N2 ) 

+ ^/(2f)Ao' - 2k B (m z (T™ ))Aa'T N2 . (15) 

If the difference between the temperatures of the met¬ 
als, Tn 1 and Tjv 2 , is not too large, then the variation 
of the damping constant Ac/ is very small, Ac/ -C a'. 
In particular, the experimental data show the following 
change of the damping constant with temperature^ if 
AT = T N1 - T N2 « 350K then |Ac/|/c/ « 0.28. For 
Tn 1 — T n 2 < 100 the relative variation of the damping 
constant is even smaller, |Ac/|/c/ < 0.1. In such a case, 
the spin pumping currents transferred into the metals al¬ 
most compensate each other, I sp 1 ss —I S p 2 - Thus, if the 
difference between the temperatures of the metals, Tni 
and T/V 2 , is not high, the IV|FIAT system becomes nearly 
symmetric, and only the spin torque current contributes 
to the total spin current. The ferromagnetic part serves 
then as a conductor which transfers the spin current from 
the hot metal to the cold one. We note, that though the 
spin pumping currents to the left and right metals do not 
contribute to the total current, I sp 1 « —I sp 2 , they lead 
to an enhanced Gilbert damping a' in Eq. (15). 

Hence, in the symmetric case, the calculation of the 
total spin current is greatly simplified, and the expression 
for the total spin current reads 


{Isz) = 2 a'k B (m z )(T N1 - T N2 ). (16) 

As follows from Eq. OSD. the spin current flowing through 
the ferromagnet mostly depends on the temperature dif¬ 
ference between the metals. However, an additional 
temperature dependence enters through the mean value 
of the magnetization (m z ), which is a function of the 
magnon temperature T™. 

In the absence of magnetic anisotropy and in the limit 
of a weak external magnetic field, M s VHq « k B Tp- : 
the expression for the spin current takes the form 


(13) 


n \ 2 T rzi Tni ~ Tn2 

{Isz) = M s VHo -—-• 


- (w® ff m z m z ) - {u)y S m z m y ). 


(17) 

















Obviously, the spin current decreases upon enhancing the 
magnon temperature. From the physical point of view 
this result is rather clear. The magnetic part conducts 
the spin torque current from the hot normal metal to the 
cold one, and the resistance of magnetic element is in¬ 
creasing with the magnon temperature T™. Therefore, 
the spin current decreases at the elevated magnon tem¬ 
perature. 

In the limit of a weak external magnetic field, 
M S VH 0 « ksT'p, the out-of-plane magnetic 
anisotropy leads to the following correction in the spin 
current formula, 


it \ 2 uj iru Tni ~ Tn2 

( Isz ) — M S VH o-—- 


2 M s VH a \ 
+ 15 h B T™ )' 


(18) 


Thus, in the case of out-of-plane magnetic anisotropy, the 
spin current always decreases with increasing the magnon 
temperature T™. This result is clear since the magnetic 
anisotropy tends to align the magnetic moment along the 
z axis, and thus increases the spin current in accordance 
with Eq. (usd, while the high magnon temperature Tp 
reduces ( m z ) and thus decreases the spin current. 

In the case of in-plane magnetic anisotropy, the expres¬ 
sion for the current reads 


(I.*) 


-a'M s VH 0 


Tni — Tn2 

rpm 

1 F 



2 M s VH a 
15 k B rp 


(19) 


Thus, the temperature dependence is similar in both 
cases. The difference concerns the role of magnetic 
anisotropy - the out-of-plane anisotropy increases the 
spin current whereas the in-plane magnetic anisotropy 
reduces the spin current. 


VI. NUMERICAL RESULTS FOR THE N\F\N 
STRUCTURE 


As in the case of the F\N structure, the analytical 
results derived above for the N\F\N system can be sup¬ 
ported by direct numerical integration of the stochas¬ 
tic LLG equation for the corresponding macrospin. In¬ 
stead of this, we go in this section beyond the macrospin 
approximation. It is well known, that the macrospin 
description breaks down for non-uniformly magnetized 
samples with the characteristic lengths exceeding several 
tens of nm’s. Beyond the macrospin formulation, the SSE 
effect can be described by introducing the local magne¬ 
tization fri [r , t ). The description can be then reduced to 
a discrete chain of magnetic moments. 

As shown in recent work of Etesami et alJ^, the spin 
current in the case of discrete chain of magnetic moments 
can be calculated by using the following recurrent rela¬ 
tions 

n 

I n = 2Aa^mf(mf_ 1 +mf +1 )-mf(mf_ 1 +mf + i), (20) 

i=i 


where A is the exchange stiffness, a is the unit cell size, 
mf, mf are the components of the individual magnetic 
moments, and I n is the site-dependent spin current. We 
note that in the case of a non-uniformly magnetized sam¬ 
ple, the spin current is not uniform along the chain. The 
chain is oriented along the x axes and the easy axes is in 
the z direction. 

The dynamics of magnetic moments fh n is described 
by coupled stochastic LLG equations, we have to solve 
numerically. For more technical details we refer to the 
recent work of Etesami et alJ^. The way how we take 
into consideration the interface effects is straightforward. 
Namely, the contact of the magnetic chain with the left 
and right metallic parts leads to modification of the 
damping constant for the first and last magnetic mo¬ 
ments, ai t N = a + 7% e ff/(47raM s ), where g e g is the ef¬ 
fective spin-mixing constant at the interfaces. Constant 
ai t N models enhancing of the Gilbert damping and is 
related to the spin pumping current. For more details 
see^i. The other point is that in the equations of mo¬ 
tion for the first and last magnetic moments, which are 
in contact with the left and right normal metals, include 
the spin torque term. 

The idea of the spin torque is that it accounts for 
the effect of the interface on the adjacent magnetic 
moment s 15 ' 25 . We assume that the ratio between the am¬ 
plitudes of the spin torque terms is proportional to the 
ratio of the temperatures of the two metals, ||/1| = 
Tl /Tr and ip = A(Tf,0,0) and /£ = A(-T*,0,0), 
where A is a phenomenological constant. 


Thus, the equation of motion for the magnetic moment 
m n can be written as 


drh n 

~dt 


7 


1 + a n 

7«n 


eff 


rh n X H n 


(1 + a n) L 

for n = 2, ...N — 1, and 
drhi t N 


m n x 


rh n x H' n 


eff 


dt 


_ 0 _ 

1 + a\ N l 

7 a i,JV 


mi,jv x Hl s N 


(1 + a\ N ) V 

7 


TOi.jv X 


TOyAT X Hl S N 


Ms a 3 
7 

Ms a 3 


myjv x 


rni ,n x ip R 


P 


mi,N x /“ft 


( 21 ) 


( 22 ) 


for the magnetic moments in direct contact with the 
metallic parts. The last two terms in Eq. E^l) stand 
for the spin currents /y n R injected from the metals to 
the ferromagnetic chain. One of them is the torque of 
Slonczewski’s type and the second one corresponds to an 
additional torque ter m 15 ' 21 ' 27 . The coefficient /? = 0.001 
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(T^—Tr)[K] 


FIG. 9: Spin current through a FM cell as a function of 
the temperature difference of the two attached metals on the 
left(L) and right (R). The red curve is in line with analytical 
result Eq. (fT51) . a' = 0.05 and T m = 1 [K], The blue curve 
corresponds to the numerical result based on Eqs. (HUmED 
and (12211 . In the FM insulator spin torque is injected from the 
both left and right metals. Spin torque injected from the left 
metal is fixed ATf; = 5 x 10 15 [fis -1 ] while spin torque injected 
from the right metal is swaped. Here A = 7.1 x 10 15 [?is _1 T' _1 ] 
is a phenomenological constant which connects temperatures 
of the metals and corresponding spin torques |/l°r| = AT£ b . 
The spin-mixing coefficient at both interfaces (left and right) 
is g e ff = 1.14 x 10 22 [m -2 ]. The temperature of the FM sys¬ 
tem is Tp — 1 [K\. The other parameters are presented in 
Tabled 


TABLE I: The parameters used for the A | FI IV system in 
numerical simulations 

DESCRIPTION VALUE 


Anisotropy constant k , [J/m 3 ] 

4.8 x 10 4 

Exchange stiffness A, [J/m] 

1.05 x 10' 

Saturation magnetization Ms, [A/m] 

1.7 x 10 6 

Gilbert damping a 

0.01 

FM cell size a, [m] 

20 x10 -9 

External field Ho, [T] 

10“ 5 

Number of realizations, R 

100 


describes the relative strength of the last torque term 
with respect to the Slonczewski’s torque. 

The results of numerical calculations are shown in the 
Fig. [9] We plotted the spin current conducted through 
the FM insulator. In numerical calculations, the FM in¬ 
sulator is modeled by a chain of coupled magnetic mo¬ 
ments. The spin current is plotted as a function of 
temperature difference between two metals T/ v — . 

As we see numerical result fits the analytical behavior 
(Eq. ijlSjU . Until the difference between metal tempera¬ 
tures is not too large, the spin current increases linearly 
with the difference T/ v — . 


VII. CONCLUSIONS 


We have studied the influence of in-plane magnetic 
anisotropy on the thermally activated spin current flow¬ 
ing through interfaces between ferromagnetic insulator 
and nonmagnetic metal. We have considered two differ¬ 
ent systems: (i) ferromagnetic insulator and nonmagnetic 
metal in a direct contact, TV IF 1 , and (ii) ferromagnetic 
insulator sandwiched between two nonmagnetic metals, 
TVlF'liV structure. In the symmetric case, K x = K y , we 
derived analytical expressions for the average spin cur¬ 
rent in several limiting situations. 

In the case of a weak anisotropy and a weak exter¬ 
nal field, M s VH 0 /k B T™ < 1 and H A M s V/2k B T™ < 1 
(this case refers to the high magnon temperature T™), 
the total spin current is given by the formula ( I sz ) = 

(Tp - T N ) (l - §f™) • The in-plane 

positive anisotropy, H A > 0, suppresses then the spin 
current, while the negative anisotropy, H A < 0, enhances 
the spin current. 


In the case of a strong magnetic field and a weak 
anisotropy, \H a M s V « k B Tp « H 0 M S V, we iden¬ 
tified a different regime, and the expression for the 

spin current reads then ( I sz ) = 2 a'k B (Tp — T A r) I 1 — 


3H a k‘B Ty 
H 0 M s VHo 


In this case, the magnetic anisotropy sup¬ 


presses the spin current, too. The situation is different 
in the asymmetric case of K x ^ K y . We find that for a 
weak magnetic field, 2Hq/ (K x + K y ) < 1, the asymme¬ 
try reduces the spin current, see Fig. [I] 


Another conclusion is that if the difference between 
temperatures TVi and T/v 2 of the two metals in the 
TV IF 1 1 TV system is not too large, the problem becomes 
effectively symmetric. The spin pumping currents flow¬ 
ing in opposite directions compensate then each other, 
Ispi = —Isp 2 , so that only the spin-torque components 
contribute to the total spin current. In this case, the fer¬ 
romagnetic element of the structure serves as a conduc¬ 
tor, which conducts the spin-torque current from the hot 
to cold metal. In the linear response regime, the spin cur¬ 
rent is linear in the temperature difference (T/vi — TV2). 
If the difference between the temperatures of two metals 
is large enough, this breaks the symmetry between the 
spin pumping currents, I sp i ^ —I sp 2 driving the system 
beyond the linear response regime. 
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Appendix A: Average values of the magnetization 
components 

The mean components of the magnetization, see 
Eq. 0, are given as 


+1 

(m z ) = ^e B I dxxe Ax ~ Bx2 I 0 (O , (Al) 
-1 

+ 1 

(1 - m 1 2 3 z ) = J dx (l — x 2 ) e Ax ~ Bx2 I 0 (Q, (A2) 

-1 


+ 1 


(m z ml) = ^e B J dxx( 1 - x 2 ) e Ax Bx2 [Io(() + h(()] 7 

(A3) 


-1 


and ( = A (l — x 2 ), while Iq(- ■ •) and /i(- • •) are the 
modified Bessel functions of the zeroth and first order, 
respectively. Apart from this, we introduced here the 
following notation: B = ^(B x +B y ) and A = ^(B x — B y ). 

In case of a weak magnetic anisotropy, B x = B y = 
B < 1, the spin current and the mean components of the 
magnetization take the form 


(1 


(I sz ) =2a'k B (T™ -T n ) 


2 \ 2L(A) ,„fdL(A) , B [A) 

"-»> = —- 4S (— + — 


)} 



2 

A 2 


(A6) 


(A7) 



+i 


(m z rriy) = ^ e B J dxx 

(l-x 2 )e Ax ~ Bx2 [Io(0 

-h( C)], 

-l 

where 

+i 

(A4) 

Z = 2we B 

J dxe Ax ~ Bx2 I 0 (0 

(AS) 


-l 


{m z {m 2 x +m 2 )) = 


4 B 


30 L(A) 
A 4 


where L(A) = coth(A) 


6 L(A) | 2 
A 2 A 

L(A) 3L 2 (A) 10 \ ( ' A9 ' ) 

~A 2 ~ + ^3 Jsj’ 

— is the Langevin function. 
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